Observation and mitigation of image distortion in high energy electron radiography* 
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Image distortion caused by the angular misalignment of quadrupole magnets in high-energy electron radiogra- 
phy has been studied systematically. We propose that the distortion originates from the coupling of the electron 
motions in the transverse directions, based on a theoretical analysis and the transfer-matrix method. The relative 
angular rotation between the second and third magnetic quadrupoles was identified as the main contributor to 
image distortion. This was verified by both a beam-dynamics simulation and experiments. Different strategies 
to mitigate this image distortion are also explored, including online magnet tuning using a higher beam energy 
and a greater magnification factor. This study provides criteria for designing experiments and paves the way for 


achieving higher image precision. 
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I. INTRODUCTION 


High-energy density physics (HEDP) is the study of matter 
under extreme pressure and temperature conditions. Under 
these conditions, the hydrodynamic response results in a high 
expansion velocity of the order of microns per nanosecond. 
This therefore requires a diagnostic method with, simultane- 
ously, a temporal resolution ranging from nanoseconds to pi- 
coseconds and a spatial resolution <1 um. Charged particle 
radiography has developed over decades [1-3]. High energy 
electron radiography (HEER), one instance of this method, 
is attracting growing interest for its combination of potential 
high spatiotemporal resolution, accessibility, and ease of ma- 
nipulation. HEER has been proposed for use in high-intensity 
heavy-ion accelerator facilities at the forefront of HEDP sci- 
ence [4, 5]. The basic idea of HEER is to use a set of mag- 
netic quadrupoles to focus the scattered electrons passing 
through an object to form a point-to-point image on the de- 
tector plane. Since its proof-of-principle demonstration at the 
Los Alamos National Laboratory [6], significant efforts have 
been made to improve its performance, e.g., by including a 
higher-brightness electron source and an advanced electron- 
imaging optics design [7, 8]. The application of high-quality 
beams with a low energy spread, low emittance, high bright- 
ness, and short bunch length generated by state-of-the-art S- 
band photoinjectors has made significant progress in pushing 
the spatial resolution to the micron level [9]. Recently, a pi- 
oneering high spatial resolution of 0.8 um was achieved with 
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a radio-frequency (rf) thermionic cathode gun at the Institute 
of Modern Physics, which provides high-flux electron beams 
with a low energy spread and ensures an excellent signal-to- 
noise ratio in the acquired images [10, 11]. 

Meanwhile, upgrading and optimizing the imaging system 
(for example by using elaborately designed quadrupoles and 
cascaded image units) also play an important role in achieving 
better spatial resolutions and expanding the scope of applica- 
tion [12, 13]. In particular, the use of a permanent magnet 
quadrupole (PMQ), with its appealing high gradient, allows a 
shorter focal length. This is beneficial for making the entire 
imaging system more compact and increasing the magnifica- 
tion factor [14-16]. Additional advantages are its lower cost, 
non-restriction from the thermal effect, and fluctuation of cur- 
rent [17, 18]. However, drawbacks of PMQs include demag- 
netization by radiation, stronger higher-order multiple fields, 
and sensitivity to misalignment, which degrade image accu- 
racy and even cause serious deterioration of the image quality 
[19-21]. For example, image skewing was observed at the 
University of California Los Angeles (UCLA) when using a 
4 MeV electron beam focused by a PMQ triplet [22]. Mitigat- 
ing these defects is therefore crucial to further improve image 
quality. 

In this paper, we report on the image distortion caused by 
the angular misalignment of PMQs in the Russian quadruplet 
(RQ) system, and propose strategies to mitigate it. Section II 
presents the experimental design and the observed phenom- 
ena. Their possible explanations are discussed and, eventu- 
ally, identified with PMQ misalignment. Section II analyzes 
the origin of the distortion using the transfer matrix and per- 
forms start-to-finish simulations to validate the analysis. We 
also evaluate the influence of different misalignment modes 
and provide an upper tolerance limit for misalignment. Sec- 
tion IV demonstrates the feasibility of online tuning to correct 
misalignments in simulations and experiments. Finally, the 
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Fig. 1. (Color online) Schematic of the experiment setup. Electrons 
penetrating the sample are focused by four PMQs placed antisym- 
metrically with their focusing direction alternating, as denoted by 
‘+A/B’, where ‘+’ denotes focusing or defocusing and ‘A/B’ rep- 
resents PMQs with different lengths. The point-to-point image is 
formed on the YAG screen and then deflected to the camera by a 
mirror. 


impacts of different factors in optical design, including the 
field gradient, beam energy, and magnification, are discussed 
and strategies are proposed accordingly. 


Il. EXPERIMENT DESIGN AND RESULTS 


The experiment was performed using the Xi’an 120 MeV 
linear accelerator. Electron beams with tunable time struc- 
tures were generated in an S-band photocathode rf-gun and 
then accelerated by two 3-m-long traveling-wave accelerat- 
ing tubes. The energy was tuned from 30 to 120 MeV. A set 
of electromagnetic quadrupoles and correctors was used to 
focus and match the beam [23]. Beams with a 400 pC bunch 
charge, 10 ps pulse duration, 1 mm.mrad normalized emit- 
tance, and 0.1% energy spread were obtained at the entrance 
of the imaging system chamber. The imaging system con- 
sisted primarily of a quadruplet, as commonly used in HEER 
for its convenience in achieving equal focal lengths and mag- 
nifications in the x- and y-directions [24, 25]. An aperture 
placed in the Fourier plane, where the scattering angle was 
transformed into a transverse distribution, blocked electrons 
with a larger scattering angle, resulting in a lower transmit- 
tance in the image plane. A 100 um-thick high-resolution 
yttrium-aluminum-garnet (YAG) screen was used to convert 
the electron profile distribution to a visible-light image. The 
image was reflected onto a lens-coupled CCD placed per- 
pendicular to the beamline using a mirror to avoid radiation 
damage. All the imaging system equipment (except the lens- 
coupled CCD) was placed in a vacuum chamber as shown in 
Fig. 1. 

The imaging optics used in our experiment were similar to 
those in [14], where four PMQs were used to form a RQ. Un- 
like the two sets of PMQs with different gradients and magnet 
lengths used in Ref. [14], our experiment used identical gra- 
dients for all the PMQs, thereby facilitating the fabrication 
process and ensuring better uniformity. As their gradients are 
fixed, the PMQ positions become the only tunable parame- 
ters. These were optimized using the high-order beam trans- 
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Fig. 2. (Color online) Beam pattern in the image plane: (a) simula- 
tion and (b) experiment. 


port code COSY INFINITY [26], which describes the fringe 
field using a six-parameter Enge function. In the experiment, 
the parameters were adjusted using several remote-controlled 
custom stages with micron-level precision. The parameters 
for a 30 MeV beam energy and magnification factors —1.85 
and —1.96 in the horizontal and vertical planes, respectively, 
are listed in Table 1. The positions were calculated from the 
center of the magnet to the object plane. 


TABLE 1. RQ parameters used in the experiment with 30 MeV 
electrons. 


Parameter Length (mm) Gradient (T/m) Position (mm) 
Al 13.47 191.4 22.4 
B1 20.15 —186.8 49.3 
B2 20.12 188.9 97.3 
A2 13.64 —189.4 124.3 


Simulations were performed before the experiments to 
validate the effectiveness of the optimized parameters. A 
Monte Carlo simulation of electrons penetrating a standard 
100-mesh copper transmission electron microscopy (TEM) 
grid (pitch width 300 um, hole width 260 um, and bar width 
40um) with a thickness of 50um was performed using 
GEANT4 [27]. The ASTRA code [28] was employed to con- 
duct a beam-tracking simulation of the imaging system using 
the measured PMQ field maps. With a 0.2 mm radius col- 
limator placed in the Fourier plane, a clear magnified image 
of the grid is expected, as shown in Fig. 2(a), 120.84 mm 
downstream of A2. This result confirms the formation of the 
point-to-point image of an object. 

However, there are always discrepancies between the ex- 
perimental and ideal designs. After fine-tuning the beam and 
the quadrupole positions, a representative single-shot image 
was obtained in the experiment, as shown in Fig. 2(b). The 
image is clear and some small contaminants are discernible. 
In spite of not applying a collimator in the experiment, the 
small magnet apertures worked equivalently to form an image 
contrast. On the other hand, the image is rotated clockwise by 
approximately 30° and the magnification factors are —2.26 
and —3.06 in the horizontal and vertical planes, respectively. 

As mentioned previously, some common factors may cause 
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Fig. 3. (Color online) Images of the target rotated by 20° (a) counter- 
clockwise and (b) clockwise. A square copper patch in the lower-left 
region is imaged as a rhombus on the screen. 


image distortion. Here, we present how these were eliminated 
in the experiment. First, because the distortion of our pattern 
appears similar to that observed at UCLA [22], the possibility 
of sample misplacement had to be excluded. Therefore, the 
sample was carefully checked before and after shooting to 
ensure it was set perpendicular to the beam longitudinal axis. 
We also tested deliberately orienting the grid at an angle of up 
to + 20° relative to the transverse plane of the beamline. The 
measurements showed that the leaning direction was indepen- 
dent of the sample rotation, thus excluding the possibility of 
image distortion induced by sample rotation, as shown in Fig. 
3. Second, the light path of the image-collecting system (in- 
cluding the mirror- and the lens-coupled CCD) was calibrated 
before conducting the experiment by imaging a standard cali- 
bration grid marker. Third, a misalignment of the beam trajec- 
tory relative to the magnetic-field center blurs and shifts the 
image. Hence, a HeNe laser copropagating with an electron 
beam was used to prealign the sample and quadruplet along 
the beamline axis. This was further mitigated by comparing 
the beam positions when the PMQs were moved into and out 
of the beamline individually. This was designed to emulate 
the transverse deflection experienced by the beam when its 
trajectory deviates from the magnetic center of the PMQs. In 
addition, optics calculations predicted unequal magnification 
along the perpendicular axes. Although the predicted effect is 
not as prominent as we observed, it causes image stretching 
instead of rotation. 


Angular misalignment arises because the quadrupole mag- 
netic field has only two symmetry planes, in contrast to the 
circular symmetrical field of a solenoid (round lens). Fur- 
thermore, the small PMQ aperture (8 mm in the present case) 


Fig. 4. (Color online) Magnetic field axis (red dashed arrows) ori- 
entation of the first magnet, with azimuthal angle 0 relative to the 
beam transverse axis (black dashed arrows). 


causes manufacturing, measurement, and alignment difficul- 
ties. Typically, the error in the magnetization direction of 
the magnet segments is approximately + 1° in manufactur- 
ing. The measurement error of the magnetic field is on the 
same scale, limited by the size of the Hall probe. Angular 
misalignment between the quadrupoles causes changes in the 
force experienced by the particle, making its trace change in 
consequence. The following section tests this hypothesis by 
investigating the influence of angular misalignment on radio- 


graphy. 


Il. THEORETICAL ANALYSIS AND SIMULATION 


A simple example of PMQs with angular misalignment is 
presented in Fig. 4. The transfer matrix of a quadrupole 
magnet with angular rotation 0 can be written 


Mọ = R-ọMo Rọ, (1) 
with 
een, 0 0 
M. = —K sin KL cos KL 0 0 
o~ inh KL |?’ 
0 0 cosh KL aeaa 
0 0 Ksinh kL cosh KL 
cos 0 0 sind 0 
iS 0 cos@ 0 sind 
9~ | -sin 0 cos@ 0 l 
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where K? = G/Bp, Bp is the magnetic rigidity related to 
the charge-mass ratio and the particle energy, G is the field 
gradient, and L is the effective magnetic field length of the 
PMQ. Since 0 < 1, Eq. 1 can be rewritten as 

Me = M+ 0AM, 


where 
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Fig. 5. Image rotation increases with the rotation angle of PMQ B1: 
(a) 0 =0°, (b) 8 = 1°, (c) 0 = 2°, (d) 0 = 5°. 


M, and My are the two upper-left and lower-right 2x2 ma- 
trices of Mo, respectively. The entire transfer matrix of the 
RQ system with the ¿th magnet rotated by 6 is obtained as 


Me! = MaA + @My_uv,AMiMy,-_o.- (3) 


Here, M, forall is the total transfer matrix without misalignment, 
Mm, —o is the transfer matrix from the object plane to the ith 
magnet, and Mr- m; is the transfer matrix from the ¿th mag- 
net to the image plane. A similar result can be obtained com- 
putationally in more complex situations. The matrix shows 
that the electron motions in the x- and y-directions are cou- 
pled, yielding nonzero upper-right and lower-left matrix ele- 
ments (in particular Mı3 and M31) and causing image rota- 
tion and blurring. 

Simulations of the rotated PMQs were performed to verify 
this prediction. For simplicity, we first considered a situation 
in which only magnet B1 was rotated by an angle 0 rang- 
ing from 1° to 5°. The results are presented in Fig. 5. The 
evolution of the image was reasonably consistent with the ex- 
perimental results. Distortion arises rapidly with an increase 
in 0, and the image is barely recognizable when @ reaches 
5°. However, an incongruity remains in so far as the mesh 
is a rhombus in the simulation but rectangular in the experi- 
mental results, as shown in Fig. 2(b). This was because the 
sample shown in Fig. 2(b) is set with a rotation that coincides 
with the direction of the image rotation, as can be verified by 
comparison with Fig. 3. 

We next discuss the upper limit of the acceptable misalign- 
ment for the purpose of guiding experimental design. The 
angular misalignment of each magnet, despite being random, 
can be represented in terms of seven modes for the RQ sys- 
tem. The direction of rotation is denoted ‘N’ (clockwise) or 
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Fig. 6. (Color online) Relation between ¢ and @ for different angular 
misalignment modes. Solid lines represent ¢ in the x-direction and 
dashed lines in the y-direction. 


‘P’ (counterclockwise), as listed in Fig. 6 and Fig. 7, where 
each letter represents a magnet from Al to A2 corresponding 
to Fig. 1(b)). The rotation angle of the image is defined as ¢ 
to quantify the rotation level: 


Yrot 


o = arctan — arctan S (4) 
x 


Trot 


where (2rot; Yrot) is the distorted position and (x, y) is the un- 
perturbed position. With the expanded expressions in Eq. 2 
(which can be derived computationally), Yrot/£rot is coupled 
in the x and y directions. This results in @ varying in different 
directions, which explains why the stripes are nonperpendic- 
ular to each other, as illustrated in Fig. 5. To compare the 
influence of different magnet rotations, the image rotation ¢ 
in the x- and y-directions of the different modes was calcu- 
lated, as depicted in Fig. 6. The most significant influence 
is displayed by the ‘PPNN’ and ‘NPNN’ modes. Setting the 
image rotation, shown in Fig. 5(b), as the upper limit, the 
acceptable angular misalignment of this RQ system becomes 
+ 0.5°, since only B1 with a 1° rotation is equivalent to the 
‘NPNN’ mode with a 0.5° rotation. A comparison of the first 
four modes shows that B1 and B2 have a stronger impact. 
This is reasonable, considering their greater lengths. Upon 
closer inspection, the image rotation is very sensitive to the 
relative rotation between B1 and B2, in all of these cases. On 
the other hand, they all show a larger ¢ in the x-direction than 
in the y-direction except for the ‘PNNN’ and ‘NNNP’ modes. 
This indicates that setting the concerned feature of the target 
along with the y-direction is recommended to minimize the 
detriment. 

To estimate the image blur quantitatively, the image error 
o is defined as 


wil (@ror— a)? + (Yro y)? 
n=1 x2 +y? 


N a 


(5) 


where N denotes the number of electrons in the beam pulse. 
The relationship between o and 0 for different modes is 
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Fig. 7. (Color online) Relation between o and @ for different angular 
misalignment modes. 
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Fig. 8. Distorted image of (a) the NPNP mode with 0 = 0.65° and 
(b) the PNNP mode with 6 = 0.95°. 


shown in Fig. 7. A clear trend shows the image error ø in- 
creasing almost linearly with 0. Setting the value of ø, as 
shown in Fig. 5(b), as the upper bound to judge the recog- 
nizability of the image, the tolerance of the entire RQ sys- 
tem becomes + 0.5°. To verify the reliability of the analysis, 
we simulated the two most complicated misaligning modes, 
‘NPNP’ and ‘PNNP’, with 0 as their upper limit, as shown in 
Fig. 8. As expected, this image is similar to that in Fig. 5(b), 
and the image on the right shows a weaker rotation, since the 
>PNNP’ mode has a lower growth in @, as apparent in Fig. 6. 


IV. STRATEGIES FOR MITIGATING THE INFLUENCE 
OF ANGULAR MISALIGNMENT 


A. Online correction 


Given the unavoidability of fabrication, assembly, and 
measurement defects, we here discuss possible remedies for 
these experiments. The above analysis identified the relative 
rotation between B1 and B2 as the most significant factor. 
This provides clues for mitigating the influence of angular 
misalignment. In the experiment, the beam pattern is moni- 
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Fig. 9. Simulated beam pattern with only B1 present in the beamline. 
The rotation angle of B1 decreases from (a) 5° to (f) 0°. 


tored instantaneously using a lens-coupled CCD. PMQs can 
be moved in and out of the beamline in stages with a preci- 
sion as high as 1 um. When only one magnet is present in the 
beamline, a fine-tuned beam with a round profile originating 
from the LINAC is stretched in the defocusing direction and 
compressed in the focusing direction of the quadrupole mag- 
netic field. A good alignment of the second magnet with the 
first causes a change in the beam pattern, but the direction of 
stretching and compression remains constant (or exchanges). 
However, an angular misalignment between the two magnets 
causes a rotation of the directions. Therefore, the magnets 
can be adjusted by observing the evolution of the beam pat- 
terns as they move into the beamline individually to mitigate 
the angular misalignment. 


This process was demonstrated by simulating changes in 
the beam pattern while only B1 remained in the beamline and 
all other parameters were set to their designed values. The re- 
sults are shown in Fig. 9. The beam pattern is oblique when 
B1 is initially rotated by 5°. After correcting for this angle, 
the beam pattern returns to the vertical orientation. This sim- 
ulation yields an angular resolution better than 1° without a 
further evaluation of the images. 


The image in Fig. 10 was obtained experimentally. Only 
a part of the beam pattern is captured owing to the limited 
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Fig. 10. (Color online) Beam pattern produced with only B1 present 
in the beamline. (a) The beam pattern leans leftward initially. (b) 
A clockwise rotation of magnet rotates the stretching direction of 
the pattern clockwise correspondingly. The bright curve on the YAG 
screen is caused by a crack. 


YAG screen size. Figure 10(a) shows the beam pattern pro- 
duced with only B1 present in the beamline and with the three 
other magnets displaced, immediately after the image in Fig- 
ure 2(b) was taken. The stretching direction leans leftward, 
confirming the existence of an angular misalignment. A man- 
ual rotation of B1 rotates the stretch direction clockwise. Ow- 
ing to the lack of remote-controlled high-precision rotational 
stages, a corrected image of the TEM grid was not obtained 
in this experiment. However, the available results suffice to 
demonstrate its feasibility. By using an image processing pro- 
gram, a higher-precision correction is achievable. However, 
that analysis and its practical application are beyond the scope 
of this study. 


B. Higher beam energy 


In addition to online tuning, the optimization of the optical 
design is crucial. The use of different beam energies enables 
the extraction of different target data, such as the mass density 
and the electromagnetic field strength [29, 30]. On the one 
hand, a wide gap between the magnets allows the convenient 
insertion of a collimator; on the other hand, a compact system 
is necessary when space is limited. Therefore, quadrupole 
spacings and lengths deserve due attention in experimental 
design. We investigated the influence of different factors on 
the angular-misalignment sensitivity of our system for a fixed 
magnification. 

According to the expanded (computationally derivable) ex- 
pressions in Eq. 3, the magnifications M4; and M33 remain 
constant when the magnet lengths and spaces are scaled in- 
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Fig. 11. (Color online) Image error and rotation plotted as functions 
of the gradient for a given beam energy and magnification. B1 was 
rotated by 1° to represent an error source. 


versely with K, which is proportional to the square root of 
the field gradient when the beam energy is fixed. Further- 
more, the error source factors M13 and M3) are related only 
to the rotation angle @ in this case. This indicates that image 
error o is independent of the gradient. The same conclusion 
for the image rotation angle @ can be derived by a similar 
analysis. We verified the prediction with a numerical illus- 
tration by setting the beam energy to 30 MeV and fixing the 
magnification at —7.65. The other parameters, including the 
lengths and distances, were recalculated and B1 was rotated 
by 1° to represent an error source. Figure 11 plots the im- 
age error o and the image rotation angle @ as functions of 
the quadrupole gradient G. The independence of these quan- 
tities to the gradient allows greater flexibility in the optical 
design because the sensitivity to angular misalignment is not 
increased by the choice of the field gradient. 

The same magnification can be achieved at different beam 
energies by tuning the distance. The beam optics are re- 
designed accordingly to satisfy the requirement of a point- 
to-point image with a fixed magnification of — 7.65, using the 
same set of PMQs with gradient of 200 T/m in Fig. 11. We 
considered quadrupole B1 with a rotation angle of 1° to rep- 
resent a distortion source. The image error ø and rotation ¢, 
plotted in Fig. 12, show a rapid drop with increasing beam 
energy. However, the trend slows down substantially for en- 
ergies beyond 200 MeV. This also explains why the distortion 
in the case of a higher beam energy of 45 MeV with the same 
set of PMQs was not as apparent in previous experiments as 
in the present results [31]. Thus, a higher beam energy is 
recommended to mitigate the influence of misalignment in 
scenarios where the penetration ability is not a major factor. 


C. Larger magnification 


Different magnifications can be achieved using a given 
magnetic lens and beam energy by tuning the object-image 
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Fig. 12. (Color online) Image error and rotation plotted as functions 
of the beam energy, for a given set of PMQs and magnification. B1 
was rotated by 1° to represent an error source. 
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Fig. 13. (Color online) Image error and rotation plotted as functions 
of the magnification, for a given set of PMQs and beam energy. B1 
is rotated by 1° to represent an error source. 


separation. This makes it convenient to locate the target and 
to zoom in to obtain a detailed characterization. Because the 
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